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I.  INTRODUCTION 

Signal  processors  used  in  radar  systems  can  efficiently  solve 
the  problems  introduced  by  the  following: 

1)  Time  of  arrival  of  the  received  signal  due  to  range 
dependency. 

2)  Large  amplitude  clutter  signals. 

3)  Incoherent  phase  of  the  returned  radar  signals. 

4)  Low  signal-to-noise  ratios  (SNR)  of  returned  signals. 

One  system  configuration  that  has  been  successful  in  these  regards 
is  the  quadrature  channel  digital  signal  processor  (Figure  1).  Prac- 
tical implementations  are  possible  due  to  reducing  the  RF  signal  to  the 
Doppler  frequency  baseband  range,  i.e. , the  video  signal.  The  range 
information  is  preserved  by  sampling  during  the  interpulse  period  at  a 
rate  which  is  determined  by  the  transmitted  pulse  duration,  i.e.,  the 
use  of  range  bins.  However,  the  target  velocity  which  is  proportional 
to  the  Doppler  frequency  is  not  preserved  by  this  type  of  processor. 
Other  processors  such  as  Doppler  filter  banks  or  fast  Fourier  trans- 
form (FFT)  systems  can  achieve  Doppler  information  at  the  cost  of  more 
hardware.  Signal-to-clutter  considerations  dictate  the  use  of  moving- 
target-indicators  (Mil),  but  this  adds  to  the  hardware  complexity. 

Since  the  processor  extracts  amplitude  Information  by  a square-law  or 
linear  uetector  while  ignoring  phase  information,  the  coherent  phase 
of  the  return  signal  Is  not  needed.  SNR  Improvement  Is  necessary  for 
improved  target  detectability  and  is  obtained  by  the  noncoherent  inte- 
gration of  the  detector  outputs. 

The  MTI  in  the  quadrature  channel  processor  must  perform  the  clutter 
canceling  (I.e.,  high  pass  filtering)  for  each  range  bln.  Two  or  three 
pulse  cancelers  are  used  in  most  radars  because  they  are  simple  to 
build;  however,  they  are  not  suitable  for  all  applications.  Other 
design  methods  have  been  devised  to  meet  various  types  of  specifications, 
e.g.,  the  original  methods  were  identified  as  delay-line  cancelers  (see 
Chapter  4 of  Reference  1).  There  are  methods  available  today  that  will 
maximize  "improvement  factor"  [2]  or  target  passband  while  meeting 
minimum  improvement  factor  specifications  [3].  The  finite  impulse 
response  (FIR)  configuration  is  must  often  used  because  of  the  desirable 
transient  response  characteristics,  i.e.,  an  N coefficient  filter  will 
reach  steady-state  output  conditions  after  processing  N sample  values. 
Infinite  impulse  response  (IIR)  filters  are  not  often  used  because  of 
their  undesirable  transient  response  properties.  This  condition  may 
change  as  a result  of  a technique  to  minimize  the  IIR  transients 
produced  by  ground  clutter  environment  through  the  use  of  an  initializa- 
tion theory  [4].  Only  FIR  filters  are  considered  here  because  of  the 
hardware  savings  that  can  be  realized. 
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There  are  two  major  alternatives  available  for  realizing  the  FIR 
filters  used  in  the  MTI.  In  radar  system  terminology  these  have  become 
known  by  names  such  as  the  moving  window  (or  sliding  window)  and  the 
fixed  window  (or  block  window)  methods.  An  N-pulse  moving  window  Mil 
will  produce  a valid  output  after  receiving  the  first  N pulses  and  sub- 
sequently one  output  for  each  additional  input.  The  fixed  window 
approach  provides  an  output  after  receiving  the  first  N pulses,  another 
output  after  receiving  the  second  set  of  N pulses,  etc.  The  fixed 
window  and  moving  window  concepts  are  depicted  for  a three  coefficient 
filter  in  Figure  2.  The  number  of  outputs  from  the  magnitude  unit 
(i.e.,  residues)  are  determined  by  the  total  number  of  pulses  available, 
Np,  and  the  number  of  pulses,  N,  used  in  une  filter.  The  moving  window 

filter  produces  residues 

- (Np  - N)  + 1 , (1) 

while  the  number  of  fixed  window  residues  are 

- Integer  [Np/N]  . (2) 

Assuming  that  Np  is  an  integer  multiple  of  N and  that  the  total  number 
I of  pulses  remains  constant  between  systems,  it  follows  that 

, ^ - »(K,„  - 1)  + 1 . (3) 

^ Two  main  considerations  are  Involved  in  the  decision  to  use  moving 

; window  or  fixed  window,  viz.,  (1)  hardware  complexity  and  (2)  SNR 

I improvement.  The  following  sections  show  that  a tradeoff  is  available 

I between  these  quantities. 


II.  HARDWARE  CONSIDERATIONS 

One  of  the  standard  methods  of  realizing  an  N coefficient 
moving  window  FIR  filter  for  one  range  bin  is  shown  in  Figure  3.  This 
system  requires  (N-1)  delay  units,  N multipliers,  and  a summing  unit 
with  N input  ports.  Alternatively,  a cascade  of  (N  - 1)  sumning  units, 
each  with  two  input  ports,  could  be  used.  The  multipliers  and  adders 
must  be  capable  of  operation  at  the  pulse  repetition  frequency  (FRF), 
i.e.,  each  time  a new  sample  is  provided  for  the  range  bin.  Typically 
this  rate  is  slow  (s:5000  Hz)  as  compared  to  the  digital  hardware  capa- 
bility (:^10  MHz).  Consequently,  it  is  possible  to  reduce  the  number  of 
multipliers  and  adders  by  time  multiplexing  single  units.  A multi- 
plexed moving  window  configuration  for  one  range  bin  is  shown  in 
Figure  4.  This  system  requires  (N-1)  delay  units  for  the  input 
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samples  and  one  delay  unit  for  accumulating  the  output  during  the  multi- 
plexing operations.  The  rate  of  ope’^ation  for  the  multiplier  and  adder 
has  now  increased  to  N times  the  PRF.  This  is  still  a very  low  rate, 
e.g.,  a nine-coefficient  filter  and  a SOOO-Hz  PRF  would  require  45,000 
multiplications  and  sums  per  second.  Consequently,  the  multiplexing 
can  be  extended  to  multiple  range  bins  as  shown  in  Figure  5.  Additional 
control  circuitry  is  required  for  routing  the  appropriate  range  bin 
samples  into  the  corresponding  string  of  (N  - 1)  delay  units.  In  addi- 
tion, the  entire  IfTI  requires  K range  bin  accumulator  units,  K • (N  - 1) 
san5>le  delay  units,  a multiplier  and  one,  and  a two-port  adder.  The 
rate  of  operation  is  now  K • N • PRF.  When  1000  range  bins  are  used 
with  a PRF  of  5000  Hz  and  a nine-coefficient  filter,  the  multiplication 

and  summing  rate  is  45  X 10^  operations  per  second.  This  exceeds  the 

easily  obtainable  rate  of  10  X 10^;  however,  special  hardware  and 
designs  are  capable  of  obtaining  these  rates. 

It  is  possible  to  simultaneously  reduce  the  number  of  delay  units 
and  the  required  rate  of  operation  of  the  multiplier/adder  by  converting 
to  the  fixed  window  approach.  Figure  6 shows  a single  range  bin  sys- 
tem; a multiple  range  bin  configuration  is  shown  in  Figure  7.  Note 
that  the  input  sanqiles  are  only  used  once;  this  eliminates  the  (N  - 1) 
delay  units  required  for  each  range  bin  in  the  moving  window  realiza- 
tion. The  range  bin  accumulators  are  still  required.  The  rate  of 
operation  is  K • PRF,  e.g.,  a 1000-range  bin,  5000-Hz  system  will 
6 

require  5 X 10  operations  per  second.  This  is  well  within  the  state- 
of-the-art  hardware  capability. 

Table  1 presents  a summary  of  the  hardware  requirements  for  the 
various  configurations  discussed  in  this  section.  The  appreciable 
savings  that  can  be  obtained  in  hardware  by  using  the  fixed  window 
configuration  is  paid  for  by  a reduction  in  the  system  performance  as 
discussed  in  the  following  section. 


III.  SNR  CONSIDERATIONS 

\ The  purpose  of  the  integrator  shown  in  Figure  1 is  to  increase 

I the  SNR  prior  to  the  threshold  comparison.  In  this  analysis,  clutter 

signals  are  considered  to  be  sufficiently  suppressed  by  the  KTI  filter 
such  that  they  are  negligible  in  comparison  to  the  wide  band  Gaussian 
noise  at  the  output  of  the  Mil  filter.  The  theoretical  analysis  is 

2 2 

performed  for  the  square-law  detector  (I  + Q ) but  has  not  been  per- 
formed for  a linear  detector  ( ) ; however,  results  obtained  by 

simulation  are  included  for  both  detectors.  It  can  be  anticipated  that 
the  two  detectors  will  have  similar  probability  of  detection  perform- 
ance characteristics  (cf..  Figure  42  of  Marcum  [5],  or  as  repeated  on 
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while  the  SNR  at  the  integrator  output  is 


SNR* 

A 


K 


sm*  . SNR* 


(8) 


For  totally  independent  residues  such  as  from  a fixed  window  MTI  it  is 
shown  (Appendix  A)  that  the  number  of  residues  integrated  is  equal  to 
the  SNR  improvement,  i.e.. 


(9) 


The  residues  are  correlated  for  a moving  window  MTI  and  depends 

on  the  total  number  of  residues  integrated,  Kj^,  and  the  coefficients 

of  the  MTI  [cf..  Equation  (A. 27)].  If  the  residues  were  totally  cor- 
related, the  Kg  would  equal  1 and  it  follows  that 

1 ^ Kg  < • (10) 


There  may  be  major  limitations  to  using  the  results  of  this  SNR 
analysis  and  extending  this  to  performance  considerations  such  as  the 
probability  of  false  alarm  and  probability  of  detection.  The  outputs 
of  the  integrator  for  a moving  window  and  fixed  window  MTI  are  described 
by  different  probability  density  functions.  Marcum  [5]  has  analyzed 
the  case  when  pulses  of  statistlcally-independent  noise  plus  signal  are 
noncoherently  integrated.  This  is  equivalent  to  analyzing  the  fixed 
window  MTI. 

An  analysis  has  not  been  performed  for  the  case  of  correlated  noise. 

An  approximation  of  the  performance  can  be  obtained  by  making  a number 

of  assumptions  and  then  using  Marcum's  results.  If  the  number  of 

residues  integrated  are  sufficiently  large,  the  central  limit  theorem 

indicates  that  the  integrator  output  will  become  Gaussian.  The  mean 

value  and  variance  of  the  output  are  then  sufficient  to  calculate 

detection  or  false  alarm  statistics.  This  suggests  that  the  Marcum 

results  for  N-pulse  noncoherent  integration  can  be  used  for  the  moving 

window  Mil  if  K is  used  instead  of  the  total  number  of  residues, 
e 

Appendix  A shows  the  derivation  of  K when  a square- law  detector 
2 2 ® 

is  used;  i.e.,  R ■ I + Q . A simulation  program  was  written  to  verify 
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these  results.  It  was  possible  to  Include  the  simulation  of  a linear 

detector,  i.e.,  R + Q^.  Comparisons  of  the  results  for  a three- 

pulse  canceler  is  shown  in  Figure  8.  It  is  noted  that  this  moving 
window  Mil  produces  approximately  three  times  as  many  residues  as  a 
fixed  window  Mil  [cf..  Equation  (3)],  but  correlation  reduces  this  to 
an  effective  increase  of  approximately  1.5  times  as  many  residues.  It 
is  also  noted  that  there  is  little  difference  in  the  linear  and  square- 
law  detectors. 

The  effect  of  the  coefficient  design  method  is  illustrated  for  the 
binomial  technique  [9]  and  the  Houts  and  Burlage  technique  [3]  in 
Figures  9,  10,  and  11.  The  relative  results  for  the  design  methods 
depend  on  the  number  of  coefficients  used  and  can  be  very  significant. 
Table  2 presents  the  mathematical  relationship  between  and  based 
on  a least-squared  error  fit  to  calculated  points. 


IV.  CONCLUSIONS 

The  implementation  of  either  a moving  window  or  a fixed 
window  MTI  in  a quadrature  processor  determines  the  amount  of  hardware 
required  and  the  SNR  in^rovement  obtained. 

The  difference  in  hardware  used  for  either  structure  is  inconse- 
quential if  only  one  range  bln  is  inq>lemented  but,  typically,  radars 
contain  a large  number  of  range  bins,  i.e.,  greater  than  250.  Due 
to  the  range  bln  requirement,  the  amount  of  hardware  for  either  struc- 
ture proves  prohibitive  unless  multiplexing  is  used.  If  multiplexing 
is  used,  the  additional  amount  of  hardware  necessary  for  moving  window 
over  fixed  window  is  not  significant  (except  for  memory);  however,  as 
previously  stated,  the  timing  restraints  for  the  multiplexed  inq>le- 
mentation  of  a moving  window  processor  is  beyond  the  present  state-of- 
the-art. 

In  the  future,  as  component  speed  increases  and  size  decreases, 
the  implementation  of  a moving  window  MTI  will  become  less  hardware 
and  speed  critical;  however,  there  will  still  be  the  problem  of  cor- 
related noise  in  the  MTI  outputs.  As  previously  presented,  the  effec- 
tive SNR  improvement  is  reduced  by  the  correlated  noise  and  is  a func- 
tion of  the  particular  MTI  used,  i.e.,  binominal,  Houts -Burlage,  etc. 

The  utilization  of  either  a linear  or  square-law  detector  affects 
the  SNR  improvement.  Even  though  the  linear  detector  is  slightly 
better,  the  difference  is  not  considered  significant.  The  linear 
detector  will  have  a dynamic  range  advantage;  Hall  [10]  has  stated  a 
collapsing  loss  disadvantage. 


NO.  OF  RESIDUES 


Figure  8.  Three-pulse  canceler. 
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Figure  11.  Nine-pulse  MTI. 


TABLE  2.  LINEAR  RELATIONSHIP  OF  K AND 
USING  LEAST  SQUARES  METHOD  ® ^ 


K 

e 

Houts-Burlage 

Binomial 

Three-pulse  Mil 

1.537  - 0.689 

(Same  as  Houts-Burlage) 

Five-pulse  Mil 

3.222  - 2.130 

1.895  - 1.128 

Six-pulse  Mil 

3.667  - 2.611 

4.721  Kp^  - 3.583 

Nine-pulse  Mil 

7.456  - 6.322 

2.469  Kp^  - 1.776 
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One  application  o£  the  analysis  is  in  specifying  a quadrature  proces- 
sing system,  l.e..  Mil  coefficients  and  number  of  transmitted  pulses. 
First,  calculate  the  SNR  improvement  for  fixed  window  and  moving  window 
systems  and  relate  this  to  the  probability  of  detection.  Then,  deter- 
mine the  cost  of  additional  hardware  for  the  moving  window  MTI.  Ulti- 
mately, a decision  must  be  made  as  to  whether  Increased  probability  of 
detection  is  justified  by  the  increased  cost  of  the  moving  window  iiq>le- 
mentatlon  over  the  fixed  window  system. 
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Appendix  A.  COMPARISON  OF  SNR  AT  OUTPUTS  OF  MOVING 
AND  FIXED  WINDOW  MOVING  TARGET  INDICATORS 


The  system  configuration  used  in  this  analysis  is  discussed  in 
Section  1 and  shown  in  Figure  1.  This  is  a baseband  system  with  quadra- 
ture channels  I and  Q.  The  kth  residue  Rj^  is  assumed  to  be  calculated 
as 


Ru  - I.,  + Q,. 


(A.l) 


The  integrator  forms  an  input  for  the  decision  element  as 


(A.  2) 


An  analysis  by  Hall  and  Ward  [7]  has  considered  the  RF  counterpart  of 
this  problem  but  does  not  use  quadrature  channels.  The  main  point  of 
their  analysis  is  to  show  the  equivalent  number  of  statistically  inde- 
pendent residues  that  are  obtained  in  a moving  window  Mil. 


1.  SNR  ANALYSIS  FOR  FIXED  WINDOW  MTI 

The  I or  Q channel  outputs  for  the  kth  residue  determination 
for  a fixed  window  MTI  are  given  by 

N-1 

y^(kNT)  - ^ C^x^(kNT- nT),{a  = I,q},  k = 1,2 , (A.3) 

n«o 

where  the  filter  has  N coefficients  C^,  ...  and  x^(  ) represents 

the  input  to  the  channel.  The  signal  component  of  the  input  and  the 
noise  component  of  the  input  are  analyzed  separately.  Since  the  output 
is  taken  at  multiples  of  N input  samples,  it  follows  that  the  output 
will  have  reached  a steady-state  value  for  the  sinusoidal  signal  com- 
ponent. The  same  statement  is  true  for  the  moving  window  Mil.  Also, 
the  signal  component  inputs  to  the  I and  Q channel  are  in  quadrature 
and  the  residue  calculation  will  yield  an  answer  that  is  proportional 
to  the  magnitude-squared  value  of  this  component.  Consequently,  only 
the  noise  component  is  considered  in  the  details  of  the  analysis. 


Assuming  vhlte  Gaussian  noise  with  variance  a and  zero  mean  value  as 
the  input  signal  in  Equation  (A. 3)  and  then  determining  the  variance 
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of  the  I or  Q channel  output,  a , gives 
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> > C C x(kNT  - nT)x(kNT  - mT) 
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(A.4) 
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where  the  statistical  independence  is  used  for  the  x(  ) terms  when 
n m. 

The  residue  term  and  the  decision  term  are  calculated  per  Equations 

(A.l)  and  (A. 2),  and  the  variance  of  A can  be  related  to  the  variance 

of  R,  as  follows 
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(A.  5) 


This  result  uses  the  statistical  independence  of  the  residues  in  the 
fixed  window  MEI.  Since  the  variance  of  A has  been  increased  by  the 
factor  of  it  follows  that  the  SNR  has  increased  by  a corresponding 

amount  since  the  signal  coiiq>onent  will  Increase  by  i.e., 


SNR 
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(A.6) 
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2.  SNR  ANALYSIS  FOR  MOVING  WINDOW  Mil 

The  analysis  for  the  moving  window  MTI  follows  the  same 
pattern  presented  in  the  preceding  section.  Two  differences  are 
encountered.  First,  the  number  of  residues  calculated  is  not  the  same 
as  for  the  fixed  window  OTI.  If  the  same  number  of  input  san^les  are 
processed  in  an  N-coef f icient  filter,  the  number  of  residues  for  the 
moving  window  MTI  are 


- »(Kj„  - X)  + 1 


(A.7) 


where  is  the  number  of  residues  calculated  by  the  fixed  window  MTI. 

The  second  difference  in  the  analysis  is  that  the  residues  are  not 
necessarily  statistically  independent.  It  is  desirable,  and  possible, 
to  perform  an  analysis  similar  to  Equation  (A. 5)  and  obtain  an  equivalent 


number  of  statistically  independent  residues,  i.e.. 


(A. 8) 


This  form  results  by  expressing  the  variance  of  A as  in  the  first  part 
of  Equation  (A. 5),  i.e.. 
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Manipulations  are  performed  to  obtain 


(A.  9) 
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(A. 10) 


where  is  the  normalized  covariance  of  the  k and  J residues.  Com- 
paring to  Equations  (A. 5)  and  (A. 8)  yields 
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(A. 11) 
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The  analysis  can  be  completed  by  a 
of  the  filter  coefficients. 


process  that  evaluates 


in  terras 


The  I or  Q channel  output  for  the  kth  residue  determination  in  a 
moving  window  MTI  is 


N-1 

yc^(hT) 

n“0 

where  the  requirements  on  h are 

h - N,  N + 1,  ....  . (A. 13) 


The  kth  residue  occurs  at  h ■ N + (K  - 1)  and  the  relationship  between 
and  is  given  by  Equation  (A. 7). 

The  determination  of  the  normalized  covariance  begins  by 
considering 


(A.  14) 


and  evaluating  the  various  terms.  First, 


(A. 15) 


Since  I and  Q are  assumed  identically  distributed,  zero  mean,  and 
statistically  independent,  then 


r2 


r2 


» 


(A.  16) 
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Using  the  procedure  for  the  expected  value  of  the  product  of  the 
squares  of  Gaussian  variables  as  suggested  on  pages  205  and  206  of 
Wozencraft  and  Jacobs  [11]  allows  a simplification  to  occur,  viz.. 
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(A. 23) 


Substituting  into  Equation  (A. 22)  yields 
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(A. 24) 


From  Equation  (A.  11)  it  is  seen  that  is  summed  over  all  k and  j. 

. . s _ ^ J J 


P|^.  will  depend  only  on  the  magnitude  of  the  differ- 

j. 

elements  and  then  recognizing  that  pj^j 
be  shown  that 


The  evaluation  of 

ence  between  k and  j.  By  constructing  a matrix  containing  pj^j  as 


p^j^  and  that  pj^j^  ■ 1,  it  can 
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The  evaluation  cintinues  by  forming 


N-1 

n ■ - C C x(kT  - nT)x(kT  - ml  + iT) 

k k+i  n m 

n»0  n^O 
N-1 

- y^  C ,C  CT^  , (A. 26 

/ ^ m-1  m 

n^O 
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where  C - 0 for  (m  - i)  < 0.  The  cr^.  term  in  Equation  (A.25)  can  be 

evaluated  as  in  Equation  (A. 4).  This  allows  the  equivalent  number  of 
residues  to  be  formulated  from  Equation  (A. 11)  as 
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(A. 27) 
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This  can  be  expressed  in  terms  of  the  number  of  residues  calculated  by 
the  fixed  window  MTl  by  using  Equation  (A. 7)  and  can  be  related  to  the 
output  SNR  by  using  Equation  (A. 8). 

A computer  program  was  written  to  compute  the  results  of  Equation 
(A. 27).  A listing  of  this  program  is  included  in  Appendix  B. 
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nnn  non  non  o n n n o n n n n n n n n 


lOGWAM  Kf  ON 


BESUVAIUBlEJm 


74/7A  0PT=1 


FTN  4.?*74355 


program  KE0N(  INRiJI«0Ul<^UT.TA>JES.TAPE6,TAPf:i ) 

THIS  PROGRA'?  CALCULATES  THF  SUP  GAIN  POSSIBLE  WHEN 
CHANGING  FOOH  rUE;)  TO  MOVING  VINOOW  MT I . THE  ANALYSIS 
ASSUTS  That  RESIOUF  = I»»?  ♦ 0®«?.  the  parameters  AR- 
KFW  = NUr'HEP  OF  RE'- lOUFS  FOP  FIXEO  WINDOW 
NC  = NUMflf R OF  FILTfr  COEFFICIENTS 
KMW  = NUHHER  OF  PESIOUFS  FOR  MOVING  WINDOW 
KE  = equivalent  NU‘'riFR  OF  RESIDUES  FOR 
MOVING  WINDOW 

DB  = KE/KFW  expressed  IN  OR 
c(  ) = coefficient  array 
CK(  ) = shifted  version  of  coef.  array 
CTC(  ) = product  of  CO  AND  CKO 
DIMENSION  COS)  .CK  OS)  .CTCOS) 

DIMENSION  .<(15)  .Y(1S) 
peal  KE 
WRITE (6.104) 

104  FORPATOHl) 

PEAD(S.99)  NC 
99  FORMAT (13) 

REAO(S.IOO)  (C(I) .I=1.NC) 

100  F0RMAT(4F20.10) 

WRITE(6.10l) 

101  F0RRAT(/-'3/..31HTHE  COEFFICIENTS  ARE  /) 
wRiTEie.io?)  (C(n.i  = i.NC) 

102  FORMAT(1X,4(F20.10.?<)) 

WR1TE(6.UI) 

111  F0R.MaT(//8X.2HNC.6X..4H  kFW.R<,3HKMW.6X.3H  KE.8X.2H0R/) 

SUM  CN*»2 

SUMCSO-0. 

DO  9 J=l,NC 
SUMCS0=SUMCS0*C ( J) »C ( J) 

9 CONTINUE 
NC:-(l=NC-l 
DO  12  KFW=2.16 

SUMMING  COVARIANCE  TERMS 

I=KFW-l 
KMWM=;NC»15.1 
KMVMiKMWM.S 
FKM=FL0AT (KMWM) /s. 

FKM=S.n»AINr (EKM) 

KMW=NC»(KFW  -1) »1 

SUMDFN=KMW 

DO  10  K=l,NCMl 

CORRELATION  PRODUCT  CM-I»CM 

SUMCTC=0. 

00  11  HP  1 = 1. NC 
CTC(MPI)=0.0 

IF( (mPI-K) .GT.OI  CTC (MP| ) =C (MOl-K) ‘CIMPl ) 

suMC  rc=suMc  rc*CTC (hp i ) 
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103 
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bestjvaiwble  copy 


OGRAM  KEON 


74/74  0PT=1 


FTN  4.?*743b'i 


CONTINUE 

RHO=SU“.r,TC/SUNCSQ 

RH0S'7=RH0»»2 

K0F.L=KMW~K 

SUMOEN=SUMDe:M*2,0»FLOAT  (KOEL ) *RHOSO 
CONTINUE 

KE=rLOAT(KMW) ••2/SUMOeN 
FRFW=FL0AT (Kr W) 
DB=10.»4LOr,10(KE/FRrw) 

WRITFIft.lOl)  NC.KFW.KMw.KE.OB 
FORMAT  (U.IR.2n0.2F10.?) 

X ( 1) =FLOAT (KMW) 

Y(n=KE 

CONTINUE 

STOP 

END 
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